How is the Diamagnetic Effect (DME) Relevant 
to Stellar Surface Phenomena? 
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Abstract. The structure of stellar atmospheres can be modified if they include embedded diverging 
magnetic fields. This is due to the presence of speed filters. Here we examine some of the relevant 
effects of speed filters, using numerical simulations of the motion of charged particles. We introduce 
the concept of "floating" particles and point out the relevance of the diamagnetic effect to the 
evolution of Maxwellian velocity distributions into non-Max wellian distributions. 
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INTRODUCTION 

The motion of charged particles in the presence of magnetic fields is a fundamental 
problem in astrophysics. Magnetic fields play significant roles in stellar atmospheres 
as vv^ell as in outer space. Diverging (geometrically) magnetic fields are common in the 
universe, so that it is essential to understand the consequences of the motion of a charged 
particle in such fields vv^hen also embedded in a gravitational field. 

For the sake of simplicity as w^ell as to obtain analytic solution, v^e mimic here the true 
magnetic field structures by using monopoles. For sufficiently small regions, a monopole 
field is a reasonable approximation to the true field (for example, field lines of a magnetic 
dipole, can be locally presented as filed lines of a magnetic monopole (see fig 1)). 
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FIGURE 1. Dipole field approximated locally by a magnetic monopole. 



SPEED FILTERS AND FLOATING PARTICLES 



We follow the motion of a charged particle, located at a height zo above a magnetic 
monopole in the presence of a gravitational field. The height zo is chosen so that the 
vertical component of the Lorentz force is equal and opposite to the gravitational force 
(XFz = 0). This choice ensures that the particle is "floating" (see fig. 2). The initial 
conditions: 

|^Wn^z|sina = mg (1) 
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\FLorentz COS a = ^ (2) 

r 

From the geometry of the problem, we obtain that 

r 

tana = — , 

zo 

and from eqs. 1 and 2, we get 

tana = — . 

gr 

Thus, 



zo = ^. (3) 
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FIGURE 2. Speed filter. 

We see from eq. (3) that at a height zo above the monopole, the "floating height" 
depends on v^. If the particle has a velocity component parallel to component to the 
field, it will oscillate around zo- For < gzo the particle falls and the diamagnetic 
effect (DME) increases until the vertical component of the Lorentz force is equal to the 



gravitational force. In this case, we have Zdown < zo- For v\ > gzo the particle rises and 
the DME decreases until the the vertical component of the Lorentz force is equal to the 
gravitational force. In that case, we have Zup > zo- We conclude that the combination of 
the DME and gravitation can operate as a speed filter. 



NUMERICAL SIMULATION 

The magnetic monopole field is 



(4) 



where /3 is a constant. 

The equation of motion of a particle (mass m, charge q) in the presence of gravitation is 
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(5) 



We follow the position of a charged particle (proton) as function of time, in the 
presence of gravitation alone as well as in the presence of both gravitation and the 
DME. The particle's initial position is chosen to be at a given height above the magnetic 
monopole. 
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FIGURE 3. Height of proton as a function of time with both dia magnetic eff ect and gravit y (red line), 
and with gravity only (blue line). Initial conditions are: Vqx = ^.5^?>Kb ■ T/rrip , Voy = Q.5^?>Kh • T /rup 
, Voz = , Zo = 3 • lO^m , and r = 5 • lO^i^ 



RESULTS 



The above results are depicted in Figs. 3 and 4, while in Fig. 5, the acceleration is also 
shown. In the presence of the DME as well as gravitation, particles fall more slowly than 
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FIGURE 4. Height of proton as a function of t ime. Initial co nditions, are the same for all graphs, except 
that of the temperature: Vqx = 0.5^?>Kb -T /rup , V^y = 0.5y^3Kb -T/mp , Vq^ = , Zq = 3 • lO^m , and 
T = 5 - 10^ K. The first frame (top left) shows height of protons as function of time in the presence of 
gravity alone. 
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FIGURE 5. Acceleration along z axis of a proton as a functio n of time wit h diamagnetic effect and 
gravity. Initial conditions are: Vqx = ^•5^?>Kh -T /rup , V^y = Q.5^?>Kb ■ T/nip , Vq^ = , Zq = 3 • lO^m , 
and r = 5 • 10^7^ 



in a gravitation field alone or can even rise (Fig. 3). We see from Fig. 5, that the effect 
of the DME increases in importance as (i.e. temperature) grows, as expected. 



CONCLUSIONS 



The simulations demonstrate that the hydrostatic structure of an atmosphere is modified 
by the DME. The values of pressure and temperature change, essentially because the 



"effective" gravity is diminished (or inverted) selectively by the faster particles of the 
normal distribution. One can expect significant changes in the ionization equilibrium. 
The restructuring of the atmosphere can cause "apparent peculiarities" in the spectra 
of stars. The relevance of the DME to stellar atmospheric structure is, thus, due to the 
modification of effective gravity, which is particle speed dependence. The essence of the 
process is the speed filter. 

In the future we will include collisions in order to follow the evolution of a 
Maxwellian distribution into a non-Maxwellian one, as required for the initial condi- 
tions in the theory discussed by Scudder (1992a, 1992b, 1994), and Anderson (1994). 

We thank S. Owocki for useful discussions and helpful suggestions. 
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